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We  investigated  the  modeling  of  counterflow  diffusion  flames  in  which  the  products  of  ammonium 
perchlorate  (AP)  combustion  were  counterflowed  against  an  ethylene  fuel  stream.  The  two-dimensional 
problem  can  be  reduced  to  a  one-dimensional  boundary  value  problem  along  the  stagnation  point  stream¬ 
line  through  the  introduction  of  a  similarity  transformation.  By  utilizing  recent  developments  in  hydro¬ 
carbon,  chlorine,  NO,  and  AP  kinetics,  we  formulated  a  detailed  transport,  finite-rate  chemistry  system 
for  the  temperature,  velocity,  and  species  mass  fractions  of  the  combined  flame  system.  A  detailed  soot 
model  is  included  which  can  predict  soot  volume  fractions  as  a  function  of  the  strain  rate  and  the  fuel 
mole  fraction.  We  compare  the  results  of  this  model  with  a  series  of  experimental  measurements  in  which 
the  temperature  was  measured  with  radiation-corrected  thermocouples  and  OH  rotational  population 
distribution;  several  important  species  were  measured  with  planar  laser-induced  fluorescence,  UV-visible 
absorption,  and  Raman  spectroscopies;  and  the  soot  volume  fraction  was  measured  with  laser-induced 


incandescence  and  visible  absorption  spectroscopy. 

Introduction 

Many  solid  rocket  propellants  are  based  on  a  com¬ 
posite  mixture  of  ammonium  perclilorate  (AP)  oxi¬ 
dizer  and  polymeric  binder  fuels.  In  these  propel¬ 
lants,  complex  three-dimensional  diffusion  flame 
structures  between  the  AP  and  binder  decomposi¬ 
tion  products,  dependent  upon  the  length  scales  of 
the  heterogeneous  mixture,  drive  the  combustion  via 
heat  transfer  back  to  the  surface.  Clianging  the  AP 
crystal  size  changes  the  bum  rate  of  such  propel¬ 
lants.  Large  AP  crystals  are  governed  by  the  cooler 
AP  self-deflagration  flame  and  bum  slowly,  while 
small  AP  crystals  are  influenced  more  by  the  hot 
diffusion  flame  with  the  binder  and  bum  faster.  Tliis 
allows  control  of  composite  propellant  ballistic  prop¬ 
erties  via  particle  size  variation. 

Previous  measurements  on  AP/binder  diffusion 
flames  in  a  planar  two-dimensional  sandwich  config¬ 
uration  have  yielded  insight  into  the  controlling 
flame  structure  [1,21,  but  there  are  several  draw¬ 
backs  that  make  comparison  with  modeling  difficult. 


First,  the  flames  are  two-dimensional  in  structure, 
making  modeling  much  more  complex  computation¬ 
ally  than  with  one-dimensional  propellant  systems, 
such  as  cyclotrimethvlene  trinitramine  (RDX)  self- 
and  laser-supported  deflagration  [3].  In  addition,  lit¬ 
tle  is  known  about  the  nature,  concentration,  and 
evolution  rates  of  the  gaseous  chemical  species  pro¬ 
duced  by  the  various  binders  as  they  decompose. 
Tliis  makes  comparison  with  models  quite  difficult. 
Alternatively,  counterflow'  flames  provide  an  excel¬ 
lent  geometric  configuration  within  which  AP/ 
binder  diffusion  flames  can  be  studied  both  experi¬ 
mentally  and  computationally. 

While  counterflow  diffusion  flames  have  been 
studied  in  recent  years  using  experimental,  theo¬ 
retical,  and  numerical  techniques  [4-10],  there  has 
been  little  work  in  which  these  tools  have  been  ap¬ 
plied  to  the  study  of  solid  propellants  in  tliis  config¬ 
uration.  While  some  preliminary  studies  of  AP  coun¬ 
terflow  flames  w'ere  made  by  Friedman  [11],  Ablow 
and  Wise  [12],  Inami  and  Wise  [13],  Wiersma  and 
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LAMINAR  DIFFUSION  FLAMES— Counterflow  Flames 


Fig.  1.  Schematic  of  experimental  configuration  for  eth¬ 
ylene  versus  AP  counterflow  flame  system. 


Wise  [14],  and  Mitani  and  Niioka  [15],  the  first  work 
tliat  combined  a  detailed  transport/finite-rate  chem¬ 
istry  model  of  AP  counterflow  diffusion  flames  with 
experimental  measurements  was  by  Tanoff  et  al. 
[16].  Tliis  study  focused  on  an  investigation  of  coun¬ 
terflow  diffusion  flames  in  which  the  products  of  AP 
combustion  were  counterflowed  against  a  methane 
fuel  stream.  By  utilizing  developments  in  hydrocar¬ 
bon,  chlorine,  NO,  and  AP  kinetics,  a  detailed  trans¬ 
port,  finite-rate  chemistry  system  for  the  tempera¬ 
ture,  velocity,  and  species  mass  fractions  of  the 
combined  flame  system  was  formulated.  Results  of 
the  model  were  compared  with  a  series  of  experi¬ 
mental  measurements. 

In  tliis  paper,  we  continue  our  investigation  of  AP 
counterflow  diffusion  flames  by  incorporating  two 
major  modifications  into  the  model  developed  in 
Ref.  [16].  First,  the  fuel  was  replaced  by  ethylene, 
and  second,  a  detailed  soot  model  was  included  in 
the  problem.  Ethylene  is  an  important  reaction  in¬ 
termediate  of  hydrocarbon  binders,  and  thus,  the 
understanding  of  tliis  system  represents  an  essential 
liierarchical  step  toward  the  detailed  analysis  of  a 
composite  propellant.  Because  composite  propellant 
formulations  are  overall  fuel-rich  and  heteroge¬ 
neous,  soot  formation  and  its  coupling  to  energy  re¬ 
lease  rates  and  radiative  feedback  to  the  surface  are 
fundamental  processes  tliat  directly  impact  propel¬ 
lant  performance  and  design.  An  accurate  soot 
model  is  therefore  a  necessary  submodel  for  any  de¬ 
tailed  composite  propellant  model. 

We  compare  the  results  of  the  computational 
model  with  a  series  of  experimental  measurements 
in  which  the  temperature  was  measured  with  radi¬ 
ation-corrected  thermocouples  and  Oil  rotational 
population  distribution,  several  important  species 


were  measured  with  planar  laser-induced  fluores¬ 
cence  (PLIF),  UV-visible  absorption,  and  Raman 
spectroscopies,  and  the  soot  volume  fraction  was 
measured  with  laser-induced  incandescence  (LII) 
and  visible  absorption  spectroscopy.  In  the  next 
three  sections,  the  experimental  and  computational 
approaches  and  the  kinetic  mechanism  are  outlined, 
after  which  the  results  are  presented  and  discussed. 


Experimental  Approach 

In  a  typical  counterflow  diffusion  flame  experi¬ 
ment,  two  laminar  plug  flow  jets  (one  fuel  and  one 
oxidizer)  are  directed  toward  each  other  and  im¬ 
pinge  in  the  middle  of  the  domain.  Properly  de¬ 
signed,  tliis  configuration  leads  to  a  one-dimensional 
flame  along  the  stagnation  point  streamline  tliat  can 
be  modeled  as  a  one-dimensional  problem.  In  ad¬ 
dition,  the  experimentalist  has  complete  control  over 
many  flame  parameters  including  fuel  chemistry, 
fuel  thermal  properties  (via  dilution,  for  example), 
and  flame  strain  rate.  In  our  case,  the  oxidizer  side 
was  not  a  jet  but  a  solid  pellet  of  AP.  Although  AP 
does  not  normally  self-deflagrate  at  pressures  below 
200  to  800  psia  (depending  on  purity),  it  was  found 
tliat  the  counterflow  diffusion  flame  was  self-sustain¬ 
ing  even  at  1  atm,  and  tliat  nearly  perfect  planar 
multiflame  structures  were  imaged  with  metliane 
and  ethylene  as  fuels.  Fuel  mixtures  can  be  selected 
to  simulate  the  decomposition  products  of  various 
binders,  or  a  single  fuel  can  be  selected  to  better 
understand  the  kinetics  involved. 

The  flame  configuration  is  illustrated  in  Fig.  1. 
The  fuel  exit  diameter  was  7.75  mm,  and  the  AP 
pellet  diameter  was  10.1  mm.  Both  the  AP  and  fuel 
flows  were  surrounded  by  nitrogen  shroud  flows 
(26.8  mm  diameter  for  the  fuel  and  23.0  mm  for  the 
AP).  The  shroud  flows  were  set  to  match  the  fuel 
and  AP  decomposition  gas  flows.  The  fuel  flow  rate 
was  controlled  with  a  calibrated  flow  meter,  but  the 
AP  decomposition  gas  flow  rate  was  determined  by 
the  AP  solid  regression  rate.  The  AP  pellet  and  the 
ethylene  jet  were  separated  by  5.0  mm,  with  an  eth¬ 
ylene  average  flow  speed  of  20.5  cm/s,  yielding  an 
AP  regression  rate  of  0.113  mm/s  and  a  correspond¬ 
ing  strain  rate  of  approximately  311  s-1.  The  AP 
surface  location  was  maintained  with  a  mechanical 
stage  tliat  offset  the  regression  rate.  The  AP  was  ul- 
traliigh  purity  (UIIP)  grade  from  Kerr  McGee  Cor¬ 
poration,  pressed  to  98.2%  theoretical  maximum 
density  (TMD).  UUP  AP  was  used  because  the 
flame  structure  and  even  the  regression  rate  are  sen¬ 
sitive  to  impurities  in  normal  research-grade  AP. 

Diagnostics  applied  to  the  counterflow  AP  diffu¬ 
sion  flame  included  radiation-corrected  thermocou¬ 
ples  and  Oil  rotational  measurements  for  the  tem¬ 
perature,  PLIF  imaging  for  OH,  CN,  and  NO 
species  profiles,  Raman  spectroscopy  for  absolute 
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major  species  profiles,  and  LII  for  the  soot  volume 
fraction.  UV-visible  absorption  spectroscopy  was 
used  to  obtain  absolute  number  densities  for  the  Oil 
and  NO,  and  for  the  soot  volume  fraction.  The  CN 
concentration  was  below  the  CN  absorption  spec¬ 
troscopy  detection  limit  of  approximately  5-10  ppm. 
The  CN  radicals  were  monitored  with  PLIF  by 
pumping  the  (0,0)  bandhead  of  the  B-X  transition 
near  388.3  nm  while  monitoring  (0,1)  emission 
around  422  nm. 

On  the  fuel  side,  where  Oil  temperature  mea¬ 
surements  were  not  feasible,  50  //m  thermocouples 
were  used  to  measure  the  temperature  distribution. 
The  AP-side  temperature  was  obtained  from  Oil 
PLIF  measurements  of  the  Rl(3),  Rl(10),  and 
Rl(14)  fines  of  the  (1,0)  A-X  OH  transition.  The 
temperature  was  then  obtained  from  a  Boltzmann 
distribution  plot. 

The  NO  profile  shape  was  obtained  via  PLIF  by 
pumping  selected  rotational  fines  (R2(7),  Q2  + 
Q12(10),  and  P2  +  P12(16)  at  236.639  nm  of  the  (0,1) 
hot  band  of  the  A-X  gamma  system),  while  monitor¬ 
ing  (0,3)  emission  around  260  nm.  The  profile  shape 
was  corrected  for  the  temperature  dependence  of 
these  ro-vibrational  fines  using  the  measured  tem¬ 
perature  profile.  The  NO  vibrational  temperatures 
obtained  from  the  absorption  spectral  fits  (to  give 
absolute  concentration  in  mole  fraction)  matched 
the  measured  OH  rotational  temperature. 

Major  species  such  as  C2II4, 02,  IICl,  CO,  II2,  and 
II20  were  measured  by  Raman  spectroscopy  using 
330  mj  of  532  nm  fight  from  a  Nd:YAG  laser  (see 
also  Ref.  [16]).  The  Raman  signal  was  focused  onto 
a  0.22  m  focal  length  spectrograph  and  optical  mul¬ 
tichannel  analyzer.  Two  spectra  were  obtained  for 
each  experiment:  one  with  vertical  laser  polarization 
and  one  with  horizontal.  All  of  the  major  Raman 
fines  monitored  disappear  when  observed  along  the 
polarization  vector,  while  interfering  signals,  such  as 
laser  induced  fluorescence  (LIF),  are  not  affected. 
Thus,  the  interference  is  removed  by  subtracting  the 
horizontal  polarization  signal  from  the  vertical.  The 
major  interference  comes  from  C2  Swan  band  LIF, 
which  does  not  occur  in  the  spectral  regions  for  the 
monitored  Raman  fines  and  occurs  only  in  limited 
spatial  regions  of  the  flame.  The  Raman  signals  were 
averaged  over  time,  between  5  and  30  s,  to  build  up 
reasonable  statistics.  (In  the  hot  regions  of  the  flame, 
the  number  density  was  so  low'  that  collected  Raman 
signals  w'ere  as  low  as  a  few  counts  per  second). 

LII  measurements  were  made  for  the  AP/ethyl- 
ene  system  to  obtain  soot  volume  fraction  and  size 
profiles.  LII  is  a  result  of  laser  heating  of  particulates 
winch  then  radiate  broadband.  The  setup  for  these 
measurements  was  the  same  as  for  PLIF,  but  the 
laser  was  tuned  off-resonance  for  molecular  transi¬ 
tions,  and  the  narrowband  optical  collection  filter 
was  replaced  with  a  broadband  filter  to  collect  more 
blackbody  incandescence.  The  particle  size  was  also 


obtained  from  the  time  dependence  of  the  LII  emis¬ 
sion  [17].  A  Gaussian-like  growth  and  decay  curve 
was  obtained  with  maximum  particle  diameters  of 
approximately  50  nm. 

The  soot  profile  wras  made  absolute  by  using  re¬ 
sults  from  line-of-sight  extinction  measurements 
[18],  A  CW  diode  pumped  Nd:YAG  laser  of  wave¬ 
length  532  nm  was  split  into  tw'o  beams,  one  of 
which  w'as  monitored  as  the  unattenuated  beam  (/„) 
and  the  other  of  which  passed  through  the  flame. 
Tliis  latter  lieam  was  attenuated  through  extinction 
by  soot  particles.  The  ratio  of  attenuated  (7)  to  un¬ 
attenuated  fight  (70),  or  the  transmission  (7/70),  w'as 
measured  with  photodiodes.  Since  the  particle  di¬ 
ameters  were  in  the  nanometer  size  range,  the  Ray¬ 
leigh  limit  [19]  w-as  valid  for  the  computations  in 
converting  measured  transmission  to  soot  volume 
fraction.  The  soot  volume  fraction  (/v)  was  obtained 
from  the  following  formula  [19]: 

/v  =  -U/KJ  ln(///0)/L  (1) 

where  L  is  the  pathlength  (in  cm)  over  which  atten¬ 
uation  took  place,  A  is  the  laser  wavelength  (cm),  and 
K a  is  the  absorption  coefficient.  In  the  Rayleigh 
limit,  K,  is  given  by  the  following  [19]: 


[(n*(l  - 


(367rn2#c) 


where  n  is  the  real  part  and  iik  is  the  imaginary  part 
of  the  index  of  refraction.  The  index  of  refraction  for 
soot  at  532  nm  was  obtained  from  Ref.  [19].  Thus, 
for  these  experiments,  /v  =  —  14.7  ln(7/7<>)/padi- 
length  in  parts  per  million. 

For  AP/ethylene  flames,  w'e  noticed  tliat/v  w'as 
lower  in  the  central  region  of  the  flame  than  on  the 
edges,  due  to  a  radial  dependence  of  the  strain  rate. 
Radially  offset  transmission  measurements  were 
made.  The  data  w'ere  too  noisy  for  an  Abel  decon¬ 
volution,  so  a  function  was  assumed  for  the  radial 
dependence  of  soot  volume  fraction,  and  tliis  w-as 
integrated  over  the  fine  of  sight  and  compared  with 
die  measurements;  dien  die  functional  parameters 
were  adjusted  for  die  best  least-squares  fit.  The  soot 
volume  fraction  distribution  obtained  had  a  mini¬ 
mum  on  die  centerline  and  increased  widi  increas¬ 
ing  radius.  The  soot  volume  fraction  value  used  to 
normalize  the  LII  profile  was  dien  an  average  of  die 
recovered  radial  distribution  over  die  soot  disk  area. 


Computational  Approach 

To  model  die  AP/ediylene  flame  system,  we  con¬ 
sider  a  laminar  reacting  flow-  stabilized  in  die  vicinity 
of  the  stagnation  plane  between  two  axisymmetric, 
counterflowing  streams.  One  stream  contains  ediyl- 
ene  and  die  odier  die  AP  decomposition  products 
as  determined  in  Ref.  [20].  The  complete  formula¬ 
tion  of  die  madiematical  model  for  solving  die  finite 
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TABLE  1 

AP  Surface  Species  (mole  fractions) 


Species 

Xk 

o2 

0.22 

n2 

0.03 

h2o 

0.38 

HC1 

0.09 

nh3 

0.08 

no2 

0.06 

n2o 

0.03 

CIO 

0.03 

Cl. 

0.03 

C102 

0.02 

HC104 

0.02 

NO 

0.01 

Fig.  2.  Comparison  between  experimentally  measured 
and  calculated  temperature  profiles  for  the  ethylene  versus 
AP  counterflow  flame. 


burner  separation  problem  with  plug  flow  boundary 
conditions  [21,22]  starts  with  the  elliptic  form,  in 
cylindrical  coordinates,  of  the  two-dimensional 
equations  describing  the  conservation  of  total  mass, 
individual  chemical  species  mass,  momentum,  and 
energy  for  the  reactive  flow  occurring  between  the 
fuel  and  the  AP  solid.  By  seeking  a  similarity  solution 
of  the  governing  equations,  we  can  reduce  the  prob¬ 
lem  to  the  solution  of  a  nonlinear,  two-point  bound¬ 
ary  value  problem  in  the  axial  direction  along  the 
stagnation  point  streamline. 

Soot  kinetics  were  modeled  as  coalescing,  solid 
carbon  spheroids  undergoing  surface  growth  in  the 
free  molecule  limit.  The  particle  mass  range  of  in¬ 
terest  was  divided  into  sections  [23],  and  an  equation 
was  written  for  each  section,  including  coalescence, 
surface  growth,  and  oxidation.  Sectional  analysis 
makes  it  possible  to  obtain  the  particle  size  distri¬ 
bution  without  a  priori  assumptions  about  the  form 


I  • 

CaH, 

0 

CO 

j  ■ 

H, 

♦ 

h2o 

► 

| _ 

Fig.  3.  Comparison  between  experimentally  measured 
and  computed  major  species  profiles  for  the  ethylene  ver¬ 
sus  AP  counterflow  flame. 


Fig.  4.  Comparison  between  experimentally  measured 
(PLIF)  and  calculated  OH  concentration  profiles  for  the 
ethylene  versus  AP  counterflow  flame.  Experimental  OH 
measurements  are  absolute  concentrations. 


of  the  distribution.  For  the  smallest  section,  an  in¬ 
ception  source  term  was  included.  The  transport 
conservation  equation  for  each  section  includes  ther¬ 
mophoresis,  an  effective  bin  diffusion  rate,  and 
source  terms  for  gas-pliase  scrubbing.  The  gas  and 
soot  equations  are  additionally  coupled  through  non- 
adiabatic  radiative  loss  in  the  optically  thin  approxi¬ 
mation.  The  inception  model  employed  here  is 
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Fig.  5.  Comparison  between  experimentally  measured 
(PLIF)  and  calculated  CN  concentration  profiles  for  the 
ethylene  versus  AP  counterflow  flame.  Experimental  CN 
measurements  have  been  normalized  to  the  peak  concen¬ 
tration  predicted  by  the  model,  but  the  spatial  positions  are 
absolute. 


Fig.  6.  Comparison  between  experimentally  measured 
(PLIF)  and  calculated  NO  concentration  profiles  for  the 
ethylene  versus  AP  counterflow  flame.  Experimental  NO 
measurements  are  absolute  concentrations. 

based  on  an  estimate  of  the  formation  rate  of  two- 
and  three-ringed  aromatic  species  (naphthalene  and 
phenanthrene)  and  is  a  function  of  local  acetylene, 
benzene,  phenyl,  and  molecular  hydrogen  concen¬ 
trations  [24].  The  contributions  from  the  inception 
processes  were  incorporated  in  the  first  sectional 
bin,  whose  lower  mass  boundary  was  set  equal  to  the 
mass  of  the  smallest  inception  species  (127  amu).  In 
the  sectional  representation  [231,  die  sectional  mass 
Ixiundaries  varied  linearly  on  a  logarithmic  scale. 
The  number  of  sections  required  for  convergence 
must  be  examined  for  each  problem  and  depends  on 
die  relative  magnitudes  of  surface  growdi  and  in¬ 
ception.  Oxidation  of  soot  by  02  and  Oil  was  treated 
as  described  in  Ref.  [25].  The  surface  growdi  rate 
was  diat  proposed  by  Colket  and  Ilall  (denoted  by 


Fig.  7.  Comparison  between  experimentally  measured 
(LII)  and  calculated  soot  volume  fraction  profiles  for  the 
ethylene  versus  AP  counterflow  flame.  Experimental  soot 
measurements  are  absolute  soot  volume  fractions. 


Fig.  8.  Computed  heat  release  rate  profile  for  the  eth¬ 
ylene  versus  AP  counterflow  flame. 


MODFW  in  Ref.  [26]).  The  mechanism  treated  both 
acetylene  addition  and  elimination,  and  hence  sim¬ 
ulates  die  rapid  fall-off  of  surface  growdi  rates  at 
elevated  temperatures.  Use  of  an  expression  widi  a 
simple  Arrhenius  form  results  in  a  gross  overpred¬ 
iction  of  soot  in  diese  flames. 

The  governing  equations  were  discretized  widi  an 
adaptive  finite  difference  algorithm  and  solved  using 
Newton’s  method  for  die  velocities,  the  tempera¬ 
ture,  die  gas-phase  species,  and  die  particle  sectional 
mass  fractions.  The  system  was  closed  widi  die  ideal 
gas  law,  and  appropriate  boundary  conditions  were 
applied  on  each  side  of  die  computational  domain. 
Local  properties  were  evaluated  using  vectorized 
and  highly  optimized  transport  and  chemistry  librar¬ 
ies  [27].  The  sectional  tliermophoretic  velocities  in 
die  free  molecule  regime  were  given  in  Ref.  [25],  as 
were  die  sectional  diffusion  velocities,  which  were 
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written  with  a  mass-weighted  mean  diffusion  coef¬ 
ficient  for  each  size  class.  In  the  optically  thin  radi¬ 
ation  model  used  in  our  calculations,  the  significant 
radiating  species,  in  addition  to  particulates,  were 
II20,  CO,  and  C02.  Given  the  length  scales  of  the 
flame  investigated,  it  is  highly  unlikely  that  self-ab¬ 
sorption  was  important.  Although  the  soot  volume 
fraction  reached  parts  per  million  levels,  the  narrow¬ 
ness  of  the  soot  shell  (1-2  mm)  mitigated  any  self¬ 
absorption  effects.  We  utilize  temperature,  species 
[20],  and  regression  rate  (velocity)  information  from 
the  experiments  for  both  the  AP  and  the  ethylene 
fuel  jet  boundary  conditions. 

Results 

In  tills  section  we  examine  the  flame  structure  of 
a  one-dimensional  AP/ethylene  counterflow  diffu¬ 
sion  flame.  Pure  ethylene  flowed  from  the  fuel  jet  at 
298  K  with  a  velocity  of  20.5  cm/s.  The  AP  surface 
was  maintained  at  its  1  atm  decomposition  tempera¬ 
ture  of  825  K.  The  species  that  emerged  from  the 
AP  surface  are  listed  in  Table  1  [20],  They  had  a  flow 
speed  of  47.8  cm/s,  which  corresponds  to  an  AP  re¬ 
gression  rate  of  0.113  mm/s  (solid  density  =  1.975 
g/cm3).  Computations  were  performed  with  a  reac¬ 
tion  mechanism  that  combined  the  mechanism  util¬ 
ized  in  Ref.  [16]  with  the  ethylene  and  benzene 
chemistry  employed  in  Ref.  [28].  The  complete 
mechanism  contains  a  total  of  105  chemical  species 
participating  in  660  reversible  chemical  reactions. 
Twenty  soot  classes  were  included  in  the  computa¬ 
tion. 

The  model  predicts  an  AP  decomposition  flame 
above  the  AP  surface,  followed  by  an  ethylene/ AP- 
products  diffusion  flame,  in  agreement  with  experi¬ 
mental  observations  of  the  OH  and  CN  profiles. 
Both  flame  structures  fie  to  the  AP  side  of  the  stag¬ 
nation  plane,  which  occurs  at  1.73  mm  above  the 
ethylene  burner.  The  stagnation  plane  is  located  on 
the  fuel-rich  side  of  the  ethylene /air  diffusion  flame. 

The  AP/ethylene  system  is  remarkably  rich  in 
structure.  Color  imaging  shows  four  distinct  flame 
regions  [29].  Near  the  AP  surface  is  an  orange  region 
containing  the  burned  gases  of  the  AP  self-deflagra- 
tion  flame,  which  has  an  extremely  short  standoff 
distance  (<50 //m).  The  second  region,  which  is  light 
blue  in  color,  is  where  Oil  radicals  first  start  to  ap¬ 
pear.  Tliis  region  is  followed  by  a  third  zone  which 
is  reddish  purple  in  color  and  contains  the  primary 
diffusion  flame  with  high  radical  concentrations.  Fi¬ 
nally,  as  one  moves  toward  the  ethylene  jet,  a  bright 
yellowish  region  appears  in  the  fuel-rich  side  of  the 
diffusion  flame,  indicating  the  presence  of  soot. 

A  comparison  between  the  computed  and  mea¬ 
sured  temperature  profiles  in  Fig.  2  illustrates  the 
overall  good  agreement  between  the  model  and  the 
experiments,  with  the  exception  of  a  slight  shift  be¬ 
tween  computed  and  measured  profiles  and  a 


steeper  slope  in  the  predicted  temperature  profile 
on  the  ethylene  side.  Both  profiles  illustrate  the  ex¬ 
tremely  rapid  rise  in  the  temperature  as  one  moves 
away  from  the  AP  surface,  with  a  slower  temperature 
decrease  as  the  ethylene  jet  is  approached.  The  rise 
in  the  temperature  as  one  moves  away  from  the  AP 
surface  is  due  primarily  to  the  high  AP  surface  tem¬ 
perature  and  the  resulting  premixed  flame  that  is 
formed  when  IIC104  and  NII3  bum  with  oxygen. 
The  temperature  rises  approximately  1700  K  within 
2  mm  from  the  AP  surface. 

Figure  3  shows  a  plot  of  comparisons  among  the 
computed  and  experimental  major  species  profiles 
of  the  AP/ethylene  system.  As  in  hydrocarbon  coun¬ 
terflow  diffusion  flames,  there  is  very  little  overlap 
in  the  fuel  and  oxygen  profiles.  The  ethylene  de¬ 
creases  monotonically  until  it  is  totally  consumed  in 
the  neighborhood  of  the  peak  temperature.  As  the 
ethylene  disappears,  significant  quantities  of  II20, 
II2,  CO,  and  C02  are  formed.  (Note  that  the  exper¬ 
imental  and  numerical  ethylene  and  water  profiles 
were  divided  by  four  and  two,  respectively,  to  fit  on 
the  same  scale  as  the  other  species  shown).  What  is 
particularly  noteworthy  is  that,  while  the  hydrogen 
and  carbon  monoxide  concentrations  are  negligible 
about  1  mm  above  the  AP  surface,  the  water  con¬ 
centration  is  non-trivial  as  one  approaches  the  AP 
surface.  This  is  attributed  to  the  water  vapor  coming 
off  the  AP  surface  and  the  water  vapor  formed  from 
the  premixed  ammonia/perchloric  acid  flame  near 
the  AP  surface  via  the  reactions  NH3  +  Oil  NII2 
+  II20  and  IICl  +  Oil  — >  II20  +  Cl.  In  addition, 
most  surface  oxychlorine  compounds  (ClOa,  HCIO4) 
disappear  within  0. 1-0.2  mm  of  the  AP  surface.  CIO, 
while  falling  rapidly  over  tliis  interval  from  its  sur¬ 
face  concentration  of  0.03  mole  fraction,  remains 
present  at  a  low  level  (0.002  mole  fraction)  for  about 
0.6-0.8  mm.  ClOII  persists  in  a  similar  fashion. 

In  Fig.  4  we  compare  the  computed  and  experi¬ 
mental  Oil  mole  fraction  profile.  The  spatial  loca¬ 
tions  of  the  two  peaks  are  almost  identical,  and  the 
peak  mole  fractions  agree  to  within  30%.  In  Fig.  5, 
a  similar  profile  is  illustrated  for  CN.  The  experi¬ 
mental  profile  peaks  within  0.1  mm  of  the  computed 
result  for  CN.  The  experimental  profile  is  moder¬ 
ately  wider  in  half-width  compared  to  the  compu¬ 
tation.  Only  relative  mole  fractions  are  reported  aue 
to  the  extremely  small  mole  fraction  peaks.  Fig.  6 
compares  computed  and  experimental  NO  mole 
fractions,  which  are  remarkable  in  their  agreement. 
The  location  of  the  peak  values  are  shifted  by  0.2 
mm,  but  the  agreement  between  the  profiles,  in¬ 
cluding  the  two-tiered  structure  with  steep  gradients 
on  the  AP  side  and  the  gradual  decay  on  the  ethylene 
side,  is  excellent. 

A  comparison  of  the  computed  and  absolute  soot 
volume  fraction  is  illustrated  in  Fig.  7.  The  figure 
indicates  that  the  soot  is  localized  in  a  very  thin  re¬ 
gion  (1  mm)  in  the  fuel-rich  portion  of  the  flame. 
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This  is  consistent  with  the  visual  observations  re- 

erted  in  Ref.  [29].  Excellent  agreement  is  obtained 
•  die  spatial  locations  of  die  computed  and  mea¬ 
sured  soot  volume  fractions,  as  well  as  die  absolute 
peak  heights  (0.41  ppm  vs.  0.43  ppm).  Relative  spa¬ 
tial  distributions  for  die  separate  processes  of  surface 
growth,  soot  inception,  and  oxidation,  as  determined 
from  die  model,  indicate  diat  die  soot  oxidation  rates 
observable  in  diis  flame  are  attributable  to  super¬ 
equilibrium  Oil  concentrations,  consistent  with  our 
previous  study  [24].  Results  not  illustrated  here  in¬ 
dicate  diat  die  peak  soot  volume  fraction  is  ex¬ 
tremely  sensitive  to  die  strain  rate  of  die  system.  An 
order  of  magnitude  reduction  in  die  peak  soot  vol¬ 
ume  fraction  can  be  obtained  by  a  moderate  increase 
of  die  strain  rate.  Tliis  is  attributed  in  part  to  (1)  die 
variation  in  peak  mole  fractions  of  key  growdi  spe¬ 
cies  such  as  acetylene  and  benzene  as  a  function  of 
die  strain  rate,  (2)  die  temperature  dependence  of 
die  inception  and  surface  growdi  processes  as  a 
function  of  die  strain  rate,  and  (3)  die  significant 
decrease  in  residence  time  as  die  strain  rate  in¬ 
creases. 

In  Fig.  8,  we  plot  die  heat  release  rate  profile 
given  by 

K 

4=2  hk<bk\Vk,  (3) 

ft- 1 

where  hk  denotes  die  endialpy  of  formation  of  die 
kdi  species,  (bk  denotes  die  molar  chemical  produc¬ 
tion  (destruction)  rate  of  die  kth  species,  and  \Vk 
denotes  die  molecular  weight  of  die  fcdi  species.  We 
note  four  distinct  features  as  one  moves  from  die 
ediylene  jet  to  die  AP  surface.  First,  at  approxi¬ 
mately  3.0  mm  from  die  AP  surface,  diere  is  an  en- 
dodiermic  region  due  to  die  pyrolysis  of  C2II4  to 
C2II2.  Tliis  region  also  corresponds  to  where  die 
soot  surface  growdi  rate  is  a  maximum.  Immediately 
diereafter,  die  heat  release  peak  has  a  local  maxi¬ 
mum  due  to  die  consumption  of  C2II2  and  die  for¬ 
mation  of  CO.  A  broad  exodiennic  peak  centered  at 
about  1.3  nun  above  die  AP  surface  results  from  die 
conversion  of  CO  to  CO2  and  II2  to  II2O.  The  over¬ 
all  maximum  heat  release  rate  occurs  just  above  die 
AP  surface  and  is  due  to  die  highly  exodiermic  de¬ 
composition  of  IIC104,  Cl03,  and  Cl02  and  forma¬ 
tion  of  IICl. 

Conclusions 

The  counterflow  geometry  has  been  demonstrated 
as  a  well-suited  configuration,  bodi  experimentally 
and  computationally,  for  studying  die  structure  of 
AP  propellant  deflagrations  at  atmospheric  pressure. 
Noil-intrusive  laser  diagnostic  techniques  have 
proven  to  be  an  accurate,  reliable  means  for  probing 
die  counterflow  flame  systems  structure,  even  over 
relatively  small  (5  mm)  separation  distances.  The 


computationally  efficient  one-dimensional  flame 
model,  with  an  appropriate  kinetics  model,  has 
proven  to  be  capable  of  predicting  all  of  the  exper¬ 
imentally  observed  features  of  die  AP  versus  ediyl- 
ene  counterflow  flame  system,  including  absolute 
concentrations  and  structures  of  NO  and  Oil,  die 
complex  spatial  activity  of  NO,  CN,  Oil,  the  soot 
volume  fraction,  and  die  temperature  profile 
diroughout  bodi  die  AP  decomposition  flame  (above 
die  AP  surface)  and  die  ediylene/ AP-products  dif¬ 
fusion  flame.  A  detailed  soot  model  reveals  a  liigli- 
sooting  region  in  die  fuel-rich  region  of  die  flame. 
Tliis  is  supported  by  experimental  LII  measure¬ 
ments.  Analysis  of  die  computed  heat  release  profile 
reveals  a  multistage  set  of  reaction  zones,  consistent 
widi  color  imaging  performed  experimentally.  The 
results  of  diis  study  show  diat  die  counterflow  ge¬ 
ometry  provides  an  excellent  configuration  in  which 
to  study  propellant  oxidizer/binder  interactions.  In 
future  work,  we  will  consider  systems  in  which  die 
decomposition  products  of  hydroxy  terminated  po¬ 
lybutadiene  are  used  as  die  fuel. 
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COMMENTS 


Hai  Wang,  University  of  Delware,  USA.  The  computed 
soot  volume  fraction  profile  seems  to  be  skewed  toward  the 
cold  fuel  jet  when  compared  to  tire  experimental  profile. 
Could  you  comment  on  whether  this  discrepancy  is  caused 
by  the  uncertainty  in  soot  transport  properties  or  the  un¬ 
certainty  in  soot  chemistry? 

Author’s  Reply.  Brownian  diffusion  of  soot  particles  is 
expected  to  be  much  less  important  than  tliermophoretic 
transport,  and  the  thermophoretic  velocities  are  thought  to 
be  predicted  accurately.  Hence,  it  is  most  likely  that  the 
discrepancy  in  the  skewness  arises  from  uncertainty  in  soot 
chemistry,  either  the  gas  phase  chemistry  affecting  soot 
growth  species,  or  the  inception  and  surface  growth  pro¬ 
cesses.  The  quality  of  agreement  shown  between  experi¬ 
ment  and  theory  was  obtained  only  when  a  surface  growth 
model  exhibiting  an  inverse  temperature  dependence  at 
elevated  temperatures  was  used,  that  is,  one  exhibiting  the 
familiar  bell-shaped  curve  in  which  the  surface  growth  rises 
to  a  maximum  in  the  vicinity  of  1800  K  and  then  declines 


at  higher  temperatures.  The  particular  surface  growth 
model  used  is  that  given  by  Collate  and  Hall  (1994),  and 
it  includes  a  process  of  elimination  of  acetylene  from  tire 
soot  radical  at  high  temperatures,  analogous  to  phenyl  rad¬ 
ical  decomposition,  and  gives  much  better  agreement  than 
a  simple  Arrhenius-type  expression  for  the  surface  growth. 


James  S.  Tien,  Case  Western  Reserve  University,  USA. 
Have  you  done  a  sensitivity  analysis  of  the  species  boundary 
conditions  on  the  AP  surface? 

Author's  Reply.  Yes,  we  have  performed  a  sensitivity 
analysis  of  the  species  boundary  conditions  for  the  AP  sur¬ 
face.  Changes  in  the  nitrogen  and  chlorine  species,  for  ex¬ 
ample,  can  affect  the  downstream  flame  structure.  As  a 
result,  we  believe  that  in  future  computations  it  will  be 
critical  for  condensed-phase  decomposition  products  to 
agree  well  with  tire  corresponding  surface  species  mea¬ 
surements  included  in  Table  1  of  our  paper. 


